AD- A 199  538 


Naval  Research  Laboratory 

Washington,  DC  20375-5000 


NRL  Memorandum  Report  6232 

Fl[r  rv.^ 


The  Effect  of  First  Reflections 
in  the  Dynamic  Tear  Test 

Virginia  M.  Gensheimer 

Mechanics  of  Materials  Branch 
Materials  Science  and  Technology  Division 

L.  M.  Brock 

Department  of  Engineering  Mechanics 
University  of  Kentucky 
Lexington,  AT  40506 

M.  I.  JOLLES 

Mechanics  of  Materials  Branch 
Materials  Science  and  Technology  Division 


September  8.  1988 


DT1C: 

ELECTE 

OCT  1  3  1988 


c-\ 


Approved  for  public  release;  distribution  unlimited 


■) 


SECURITY  CLASSIFICATION  of  tH1S  PAGE  _ 


REPORT  DOCUMENTATION  PAGE 


form  Approved 
OMB  No  Q7G4-0188 


la  REPORT  SECURITY  CLASSIFICATION 

UNCLASSIFIED 

1b  RESTRICTIVE  MARKINGS 

2a.  SECURITY  CLASSIFICATION  AUTHORITY 

3  DISTRIBUTION  /  AVAILABILITY  OF  REPORT 

Approved  for  public  release;  distribution 
unlimited . 

2b.  DECLASSIFICATION /OOWNGRAOING  SCHEDULE 

4  PERFORMING  ORGANIZATION  REPORT  NUMBER(S) 

NRL  Memorandum  Report  6232 

5  MONITORING  ORGANIZATION  REPORT  NUMBER(S) 

6a.  NAME  OF  PERFORMING  ORGANIZATION 

Naval  Research  Laboratory 

6b  OFFICE  SYMBOL 
(if  applicable) 

Code  6382 

7a  NAME  OF  MONITORING  ORGANIZATION 

6c  ADORESS  (C  .-y,  State,  and  ZIP  Code) 

Washington,  DC  20375-5000 

7b  ADDRESS  (C/fy,  State,  and  HP  Code) 

8a.  NAME  OF  FUNDING /SPONSORING 
ORGANIZATION 

Office  of  Naval  Research 

8b  OFFICE  SYMBOL 
(If  applicable) 

0NR 

9  PROCUREMENT  INSTRUMENT  OENTif  iCATiON  NUMBER 

Arlington,  VA  22217 


PROGRAM 
ELEMENT  NO 

61153N22 


PROJECT 

TASK 

NO  RR022- 

NO 

o 

1 

OO 

WC"K  .NIT 
ACCESSION  NO 

DN480-509 


1 1  TITLE  (Include  Security  Classification ) 

The  Effect  of  First  Reflections  in  the  Dvnamic  Tear  Test 


12.  PERSONAL  AUTHOR(S) 

Gensheimer,  V.M. ,'  Brock,*  L.M.  and  Jolles,  M.I. 


13a.  TYPE  OP  REPORT 


13b  TIME  COVERED 
PROM _  TO 


14  DATE  OP  REPORT  (Year.  Month.  Oay) 

1988  SeDtember  8 


16  SUPPLEMENTARY  NOTATION 

*Department  of  Engineering  Mechanics,  University  of  Kentucky,  Lexington,  KY  10506 


17 

COSATI  COOES 

FIELD 

GROUP 

SUB-GROUP 

18.  SUBJECT  TERMS  (Continue  on  reverse  if  necessary  and  identify  Oy  block  number I 

Dynamic  fracture  Dynamic  tear  test 

Finite  element  Elastic  response 

Stress  waves  CA:  r-  r  -  ~V  — 


19  ABSTRACT  (Continue  on  reverse  if  necessary  and  identify  by  block  number) 

"  The  dynamic  tear  (DT)  test  if  often  used  to  study  dynamic  fracture.  Previous  work  has  indicated  that 
the  first  impact  waves  arriving  at  the  specimen  notch  may  not  create  a  stress  state  severe  enough  for  fracture. 
In  this  article,  a  finite  element  program  for  the  exact  dynamic  simulation  of  an  actual  DT  test  is  first  com¬ 
pared  with  an  analytical  solution  for  the  time  interval  prior  to  reflected  wave  appearance.  The  finite  element 
results  are  then  used  to  study  the  exact  stress  field  near  the  specimen  notch.  This  study  demonstrates  that  is 
is  the  arrival  of  the  first  reflected  waves  at  the  notch  which  generate  the  stress  field  necessary  for  fracture. 
These  results  are  then  used  to  estimate  a  critical  dynamic  stress  intensity  factor  value  for  the  DT  test.  /  , 


20  DISTRIBUTION  /  AVA1LA8ILITY  OP  ABSTRACT 

ffluNCLASSIFlED/UNL.MITED  □  SAME  AS  8PT  Q  OTIC  USERS 


22 a  NAME  Of  RESPONSIBLE  INDIVIDUAL 

Virginia  M.  Genshe 


DO  Form  1473,  JUN  36  Previous  editions  are  oo s^<e(* 

S/N  0 1 02 -LF-0 14-6603 


21  ABSTRACT  SECURITY  CLASSif 'CATION 

□  one  users  UNCLASSIFIED  _ 


22b  TELEPHONE  (include  Area  Code)  22 c  'yV9P_ 

(202)  767-9027  Code  n382 


ous  editions  are  oos«^eie  _ 'EC..Pi7y  C..asS<F'CA  r  ON  > 


l 


CONTENTS 


INTRODUCTION  .  1 

REFLECTION-FREE  SOLUTIONS  .  3 

EXPERIMENTAL  IMPACT  LOAD  .  7 

FINITE  ELEMENT  PROGRAM  OUTLINE  .  8 

FINITE  ELEMENT  CALCULATIONS  .  9 

COMPARISON  OF  ANALYTICAL  AND  FE  SOLUTION  .  10 

ESTIMATE  OF  CRITICAL  STFESS  INTENSITY  FACTOR  .  11 

BRIEF  DISCUSSION  .  11 

ACKNOWLEDGEMENT  .  12 

REFERENCES  .  13 


THE  EFFECT  OF  FIRST  REFLECTIONS 
IN  THE  DYNAMIC  TEAR  TEST 


Introduction 

Dynamic  brittle  fracture  studies  are  often  (Timoshenko,  1913;  Mason,  1936; 
Lee,  1940;  Eringen,  1953;  Nash  and  Lange,  1969;  Zehnder  and  Rosakis,  1986)  based 
on  the  dynamic  tear  (DT),  or  three-point  bend,  test.  In  this  test,  a  mass  im¬ 
pacts  the  unnotched  side  of  a  simply-supported,  notched  rectangular  bar,  and  the 
resulting  elastic  wave  diffraction  at  the  notch  causes  fracture  there.  As  in¬ 
dicated  schematically  in  Fig.  1,  test  dimensions  can  be  chosen  so  that  the  mass 
strikes  more  or  less  uniformly  across  the  bar  width,  and  complete  symmetry  exists 
with  respect  to  the  notch  plane.  8ecause  the  notch  gap  is  small,  the  wave- 
induced  stress  state  at  the  notch  is,  then,  approximately  that  of  a  Mode  I  crack. 

The  elastic  waves  are  not  treated  explicitly  in  studies  such  as  those 
referenced  above,  perhaps  in  part  because  of  the  mathematical  difficulty  inherent 
in  problems  whose  solution  requires  the  coupling  of  dilatational  and  rotational 
waves  in  geometries  which  feature  corners  and  characteristic  lengths,  see  for 
example  Miklowitz  (1982).  However,  as  pointed  out  by  Kolsky  (1973),  and  demon¬ 
strated  analytically  by  Brock  and  Rossmanith  (1985)  and  Brock  (1986),  the 
arrival  of  specific  waves  can  have  a  pronounced  effect  on  a  crack  tip  stress 
field.  As  a  first  step  in  accounting  for  the  elastic  waves  arising  in  the  DT 
test,  a  transient,  two-dimensional  analysis  was  performed  by  Brock,  Jo 1 1 es  and 
Schroedl  (1985),  hereafter  known  as  BJS.  The  analysis  used  an  impact  time 
history  determined  experimentally  by  Nash  and  Lange  (1969),  hereafter  known  as 
NL,  and  yielded  exact  solutions  for  times  prior  to  the  arrival  of  reflected 
waves.  The  bar  cross-section  at  the  notch  was  found  to  go  into  compression 
wnen  the  first  dilatational  wave  from  the  mass  impact  arrived.  The  normal 
stress  on  the  cross  section  then  peaked  in  compression,  and  subsequently 
became  tensile.  This  tensile  state  occurred  just  prior  to  the  arrival  of  the 
first  impact-generated  rotational  wave,  and  grew  rapidly  thereafter.  However, 
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che  occurence  coincided  roughly  with  what  would  have  been  che  arrivals  or  me 
first  reflected  waves  had  they  been  included  in  the  analysis.  Moreover,  ever, 
without  them,  che  Mode  I  dynamic  stress  intensity  factor  did  net  reach  any 
value  thought  (Shockey,  Kalthoff,  Klemm  and  Winkler,  1983)  to  be  critical  for 
brittle  fracture  until  well  after  the  experimentally-decermined  (NL)  fracture 
initiation  time. 

These  results  indicated  that  reflected  waves  probably  generate  the  stress 
state  necessary  to  cause  brittle  fracture  at  the  notch,  and  that  the  arrival 
of  the  earlier  reflected  waves  is  an  important  event  in  the  critical  stress 
state  generation. 

In  order  to  move  beyond  the  conjectural  nature  of  these  results,  and  to 
shed  light  on  the  critical  value  of  the  dynamic  stress  intensity  factor,  this 
article  extends  the  work  of  BJS  by  performing  an  analysis  of  the  ML  DT  test 
which  can  account  for  the  effect  of  both  impact  and  reflected  waves.  As  in 
BJS,  the  test  is  modeled  as  a  rectangular,  simply  supported  bar  subjected  to 
a  uniform  line  load  across  the  bar  width.  It  is  assumed  that  a  £lane  strain 
state  exist,  that  the  bar  remains  fixed  to  che  supports  and  that  the  notch  is 
essentially  a  Mode  I  crack  tip.  Once  again,  che  actual  impact  history  obtain 
in  NL  is  used  to  determine  the  line  load. 

A  finite  elemenc  (FE)  program  is  used  to  solve  the  full  transient  proble 
and  results  are  compared  with  those  obtained  by  an  analysis  similar  to  that 
used  in  BJS  in  the  time  interval  prior  to  the  arrival  of  reflected  waves; 
that  interval  the  anlaysis  is,  of  course,  exact.  'While  the  FE  program  does  no 
treat  individual  waves,  the  mesh  and  code  adopted  satisfy  rather  strict  error 
bounds,  and  guarantee  accurate,  numerically  stable  results. 

From  the  FE  solucion,  stress  field  histories  are  obtained,  and  are  com¬ 
pared  with  the  reflection- free  histories  to  confirm  the  observat-ons  if  BJS 
concerning  wave  reflection  effects.  An  FE-based  dynamic  stress  intensity 


factor  is  tnen  used  to  estimate  the  value  critical  for  fracture. 

In  the  next  section,  the  reflection-free  case  is  described,  and  expres¬ 
sions  for  the  normal  stress  components  presented.  In  8JS,  only  the  dynamic 
stress  intensity  factor  was  actually  given.  In  the  sequel,  the  particular 
FE  process  used  for  the  complete  problem  is  described,  and  some  results 
presented.  The  FE/ reflection- free  comparison  and  the  estimation  of  the 
critical  dynamic  stress  intensity  factor  follow. 

Reflection-free  solutions 

Consider  the  cracked  half-plane  shown  in  Fig.  2:  The  assumptions  that  plan 
strain  exists  in  the  bar  of  Fig.  1  and  that  the  mass  impact  can  be  reoresented 
by  a  time-varying  line  load  P  imply  that,  until  the  elastic  impact  waves  are 
reflected  from  the  ends  and  notched  side  of  the  bar,  the  total  bar  disturbance 
can  be  represented  in  the  geometry  of  Fig.  2.  In  Fig.  2,  we  see  that  the 
impact  waves  are  those  which  radiate  from  the  impact  site,  and  those  which 
result  from  the  subsequent  diffraction  at  the  crack  tip.  There,  h  represents 
the  ligament  length  at  the  notch  plane  in  the  bar.  The  parameter  s  is  the  time 
t  after  impact  initiation  multiplied  by  the  dilatational  wave  speed,  while  m 
is  that  speed  non-dimensionalized  with  respect  to  the  rotational  wave  speed. 

The  half-plane  is,  like  the  bar,  assumed  to  be  linearly  elastic,  isotropic 
and  homogeneous. 

If  we  now  further  restrict  our  time  of  interest  to  exclude  the  reflection 
oc  the  -impact  waves  from  the  impact,  i.e.  half-plane,  surface  itself,  then 
the  solution  at  any  point  will  be  the  superposition  of  the  two  elastic  impact 
wave  ensembles  shown  in  Fig,  2. 

In  essence,  this  solution  was  obtained  by  SJS  by  adapting  a  scheme  intro¬ 
duced  by  Freund  (1974)  for  dynamic  crack  problems  which  exhibit  a  characteris¬ 
tic  lenoth  in  the  crack  plane.  However,  as  noted  earlier,  only  the  dynamic 


stress  intensity  factor  was  actually  given,  and  to  calibrate  the  FE  study, 
complece  expressions  for  the  stresses  are  required.  These  are  presented, 
therefore,  below: 

The  impact  wave  ensemble  shown  in  Fig.  2a  constitutes  the  solution  to  a 
two-dimensional  Lamb's  problem  (Lamb,  1904)  for  a  time-varying  line  load 
(point  force  in  the  plane).  This  problem  has  been  considered  both  in  BJS  and 
by  Achenbach  (1973).  Extension  of  those  results  and  rewriting  them  in  terms 
of  the  more  convenient  coordinate  systems  indicated  in  Fig.  2a  yields  for 
x>0,y<h  the  normal  stress  components 


uere  Re(  )  and  Im(  )  denote  the  real  and  imaginary  parts  of  a  complex  function, 
(*)  denotes  s-differentiation,  and  a  subscript  j  implies  that  the  terms  in 
parentheses  are  functions  of  (p.,q.),  where 

J  J 

d2q]  =  iTx  +  (h-y)pr  =  /(r2-d2)  (3) 


d  -  /tx2  +  (h-y)2] 
while 

2  2  2  2  i 

d  q2  =  “  irx  +  (h*y)P2»  ^2  “  d  1 5  p  ~  md 


(4) 

(5) 


for  d>mx.  For  mx>d,  nowever, 

p  =  x  +■  (h-y ) /(m2-l )  >mr  (6) 

which  requires  that 
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i  d  ~  q  2  =  xx  +  (y-h)p2 

when  md>x.  In  eqns  (1)  and  (2), 


a  =  /( 1+q2) ,  b  =  /(m2+q2) 


(8) 


T  =  b2+q2 , 


R  =  4q2ab-T2 


(9) 


where  R  is  the  Rayleigh  function,  and  it  is  understood  that  Re(a,b)>0  in  the 

q-plane  with  branch  cuts  Re(q)=0  and  jlm(q)|>l,  jlm(q)|>m,  respectively. 

The  impact  wave  ensemble  in  Fig.  2b  shown  radiating  from  the  crack  tip  is 

due  to  the  diffraction  of  the  ensemble  considered  above.  In  8JS,  it  was 

found  that  tne  former  could  be  obtained  as  a  weighted  superposition  of  the 

latter.  Extension  of  these  results  yields,  therefore,  expressions  for  the 

normal  stress  j  when  x>0  in  the  form 
x 


=  ./!  r*lT'r) 

r+h  ‘l 


s 


-/ 

mn+r 


jF  (T“r) 


s 


/ 

mh+p 


r-(x-p) 

fn  F  •  QS-,dkP(s-T)dT 
1  c  ' 

r(x-p) 

Jn  F2  •  QS2dkP(s-r)dT 


(10) 


where  F .  =  F .  (p  .  ,q  .)  and 
J  J  J  J 


r 


2 


iy(x-kh)+xpi , 


P-|  =  «/[(x-kh)2-r2] 


(ID 


r  =  /(x2+y2) 
wni  le 

2  2  2  2 
r  d2  =  iy(t-kn)+xp2>  p2  =  /|(x-kh)  -m  r  | ,  c  =  mr 


(12) 


,13) 


for  r>-my.  "or  my>-r,  however, 
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2 

o  =  -y+x/(m  - 1  )  >mr 

(14) 

wnich  requires  that 

i r“q2  =  -y(r-kh)-xp2 

(15) 

when  mr>T.  The  functions  Fj(p,q)  themselves  are  defined  as 

2(’-m2)pF1  =  -Re(T2Mj 

(16) 

2('-m2)pF2  =  Re(4q2b2M+) 

(17) 

where 

-  _  /0-iq) 

"+  (n-iq)(k-iq)G+ 

(18) 

.  1  y  du  ^4u2ctB> 

‘  '  J,  u+Tq  ta  (~p_) 

(19) 

a  =  /( u2-l ) ,  8  »  /(m2-u2) ,  K  =  m2-2u2 

(20) 

and  IT.a.b)  are  given  again  by  eqns  (8)  ana  (3),  with  (a,b)  defined  in  the 
same  cut  q-plane„  We  can  obtain  the  expression  for  o  from  eqn.  (10)  by 
replacing  F-,  with  its  negative,  and  by  replacing  eqn.  (16)  with  the  definition 

2(i-m2)pF1  =  -Re[(m2-2a2)TM+]  (21) 


:or  both  a  and  a,  ,  the  functions  Q(k),  S.(k)  are  given  by 
x  y  j 


QO 


/(1+k) 
(k+n )G(k ' 


1  rm  du 


inG(k)  ..if  JH- tan*1 
1  u  * 


,23! 


/( k  -1)S1  * 


kMK 

4k2ct/(k2-l+m2)-M2 


M  =  2k2-2+m2 


(24) 
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3  =  /(k2-m2) 


S 


2 


?  ?  ? 

46k“» ( k  -mw+l ) 

- 5 - o - J  T  > 

4B^k/(kc-m-+l)-r 


;  25 ; 


Finally,  1/n  is  the  Rayleigh  wave  speed  non-dimensional i zed  with  respect  to 
the  dilatational  wave  speed;  it  is  obtained  from  the  zeroes  q  =  tin  of  R. 


Experimental  impact  load 

Impact  tests  leading  to  the  brittle  fracture  of  A517  steel  bars  in  the 
Fig„  1  configuration  have  been  performed  (NL),  and  the  time  history  of  the 
impact  force  measured.  A  plot  of  this  history  is  given  in  Fig.  3.  In  BJS, 
this  time  history  was  found  to  be  accurately  represented  by  a  line  load  P, 
uniform  across  the  bar  width,  with  the  time  response 

?(s)  1261 

s  =  0.35m,  c,  «  5790  m/sec,  Fn  =  1 .4( 103)MN/m-sec  (27) 

where  s=0  is. the  start  of  impact,  sp  corresponds  to  the  instant  of  peak  impact 
force,  end  c-j  is  the  dilatational  wave  speed  in  the  bar.  In  addition,  the 
dimensionless  wave  speed  parameters  for  the  bar  were  found  to  have  the  values 

m=1.34,  n  =  1.987  (28) 

For  the  NL  DT  test  specimen,  h=28.6  mn,  so  that  sp-12h,  and  the  time  of  fracture 
is  s**14h. 

Substitution  of  eqn.  (25)  into  the  expressions  given  in  the  previous 
section  yields,  therefore,  the  normal  stresses  at  points  within  the  bar  wnich 
are  free  of  reflected  waves.  These  expressions  are  exact  solutions  of  the 
governing  equations,  yet  involve  the  actual  NL  impact  time  history.  It  should 
be  noted  that,  in  terms  of  h,  eqn.  (26)  is  valid  for  0<s<14h,  the  first  impact 
waves  arising  from  notch  diffraction  reach  the  impact  region  at  s  = 2h  and, 
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for  she  oars  used  in  8L,  the  first  reflected  waves  to  reach  the  notch  are  from 
notched  side  of  the  N|_  bar  and  arrive  at  s~1.39h.  Therefore,  eqn.  (25) 
itself  includes  the  effects  of  wave  reflections  and,  strictly  speaking,  the 
expressions  for  and  o^  are  valid  for  0<s<1.89h. 

Finite  Element  Program  Outline 

The  FE  code  ABAQUS  (Hibbit,  Karlsson  and  Sorensen,  1984)  is  used  to 
computationally  simulate  the  DT  test.  Symmetry  with  respect  to  the  notch  Diane 
requires  that  only  half  of  the  specimen  by  modeled.  This  half  is  represented 
by  a  mesh  consisting  of  1539  nodes  and  1456  CPE4H  elements.  The  latter  are 
identical  plane  strain  square  elements,  1.588  mm  on  a  side,  with  linear  dis¬ 
placement  interpolation,  and  an  independent  hydrostatic  stress  component 
(Hibbit,  Karlsson  and  Sorensen,  1984b).  This  component  is  coupled  to  the 
constitutive  relationship  for  the  specimen  material  by  LaGrange  multipliers. 

The  use  of  such  a  hybrid  element  is  required  in  regions  of  large  deformation 
in  order  to  eliminate  artificial  stiffening  caused  by  the  propagation  of 
unrealistic  nodal  displacement  constraints  through  the  mesh.  In  order  to 
avoid  localized  element  crushing,  besides,  the  line  load  defined  by  eas.  (26) 
ana  (27)  is  actually  represented  as  a  uniform  pressure  spanning  two  elements. 

It  should  be  noted,  further,  that  the  mesh  does  not  include  special 
elements  which  could  model  a  crack  tip,  i.e.  notch,  singularity. 

In  reoard  to  its  execution,  ABAQUS  uses  an  explicit  integration  scheme, 
whereby  quantities  at  time  t.  are  determined  as  functions  of  the  same  quantities 
at  time  t. ^  .  To  eliminate  numerical  instabilities  associated  with  time-step 
chanaes  within  the  integration  scheme,  the  Hi Iber-Hughes-Taylor  single-parameter 
operator  is  employed,  with  a  damping  factor  of  -0.05. 
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In  order  to  avoid  unwieldly  trial-and-error  convergence  checks  to 
determine  the  permissable  time-step  size,  ABAQUS  provides  an  analytical  method 
The  solution  is  calculated  at  time  t  •  +  ^  based  on  an  initial  time-step  size. 

The  residual  of  the  momentum  balance  equations  are  then  determined  at  time 
t.  +  At/2.  The  relative  value  of  this  half-step  residual  is  generally  in¬ 
versely  proportional  to  solution  accuracy,  so  that  the  original  At  can  be 
adjusted  to  control  solution  stability.  The  error  bounds  for  the  half-step 
residuals  are  based  on  the  expected  nodal  forces.  The  original  time-step 
is  reduced  and  the  solution  at  t^  recalculated  if  the  half-step  residuals 
are  not  within  the  prescribed  limits. 

The  maximum  possible  time-step  size  At  is  based  on  elemental  characteris¬ 
tic  dimensions  and  material  wave  speed.  In  the  present  analysis,  the  initial 
time-step  is  chosen  to  be  the  maximum  possible  time  step. 

At  =  d*/c2  (29) 

where  c2  is  the  rotational  wave  speed  and  d*  is  the  minimum  dimension  of  the 
smallest  mesh  element.  This  time-step  is  consistent  with  that  employed  in 
another  time-dependent  three-point  bend  analysis  (Nakamura,  Shih  and  Freund, 
1986).  No  reduction  in  time-step  was  found  to  be  necessary  for  solution 
stability  in  the  present  analysis. 

Geometric  nonlinearities  due  to  large  strains  and  rotations  are  included 

in  the  present  analysis.  However,  material  nonlinearity  is  not  required,  as 

the  properties  for  A517  steel  listed  here  indicate: 

Young's  modulus  E:  206.8  GPa 
Shear  modulus  y:  75.8  GPa 
Poisson's  ratio  v:  0.27 
dilutational  wave  speed  c^ :  5790  m/sec 

rotational  wave  speed  c2:  31*17  m/sec 

Finite  Element  Calculations 

As  the  list  of  properties,  Fiq,  3  and  eq.  (27)  show,  the  Nl-based  line 
load  reaches  its  peak  at  t=6Q.46  ysec(s3s^)  where,  of  course,  t=G(s=0)  is 
the  instant  of  initial  impact.  For  purposes  of  incorporating  the  effects  of 
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at  least  the  first  dilatational  and  rotational  waves  and  notch  surface  Rayleian 
waves  (Brock  and  Rossmanith,  1985)  generated  by  reflections  from  the  notcneo 
surface  of  the  specimen,  it  is  sufficient  to  model  the  DT  test  from  t=u  to 
t=28.5  ysec.  As  examples  of  the  FE  results  generated,  the  four  non-zero 
stresses  along  the  line  y=0,x>0  from  the  crack  tip  near  the  arrival  times  of  the 
first  impact  wave  (5.32  usee)  and  the  first  reflection  (9.88  usee)  are  shown  in 
Figs.  4  and  5,  respectively.  The  variations  with  t  of  the  Mode  I  openinq 
stress  a  small  distance  (1.32  ran)  from  the  crack  tip  along  the  radial  lines 
0=0  and  d-n/2  are  shown  in  Fig.  6. 

To  study  the  Mode  I  dynamic  stress  intensity  factor  Kj,  the  classical 
(Sneddon  and  Lowengrub,  1969)  result  for  a  stationary  semi -i nfi ni te  crack  is 
used: 

Kj  =  1  im/(2irr)  ^  ■yy  ,  f(6)  =  cos  tj^-j-q)  [l+sir>^(-j’-0)sin'j(-j-9) ]  (30a, b) 

r-*0 

Because,  the.  FE  mesh  does  not  contain  a  singularity-producing  element,  Kj 
cannot  be  obtained  directly  from  the  FE  data.  However,  good  agreement  between 
FE  and  analytical  results  has  been  shown  when  the  limiting  process  of  eq.  (30a) 
is  performed  by  graphical  extrapolation  (Broek,  1982).  From  such  extrapolations 
at  various  times  during  the  interval  0<t<28.5  usee,  it  is  possible  to  develop  the 
Kj  time  history  shown  in  Fig.  7. 

Comparison  of  Analytical  and  FE  Solution 

A  partial  calibration  of  the  FE  program  can  be  obtained  by  comparing  its 
numerical  results  with  those  qenerated  by  the  analytical  solution  given  above. 

A  typical  and  rather  severe  comparison  is  shown  in  Fig.  8,  where  the  crack- 
opening  stress  on  the  notch  plane  (  =  /2)  very  near  (r=1.32  m)  the  crack 
tip  is  plotted  vs.  time.  There  is  no  disturbance  at  this  point  until  the  first 
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impact  wave  arrives  at  t=4.9  psec.  The  first  reflected  wave  arrives  at 
t=9.4  usee,  at  which  time  the  analytical  results  are,  of  course,  no  longer 
exact.  In  view  of  this.  Fig.  8  shows  excellent  agreement  until  this  time 
and,  indeed,  somewhat  beyond  it.  In  particular,  both  solutions  show  clearly 
the  intial  compression  of  the  notch  region.  Fig.  8  also  shows  that  after 
t— 1 1  psec,  the  analytical  results  begin  to  seriously  underestimate  the  magni¬ 
tude  of  the  notch  stress  state.  This  clearly  indicates  the  importance  of  the 
reflected  wave  arrivals  in  determining  this  state  and,  in  turn,  confirms  the 
conjectures  of  BJS. 

Estimate  of  Critical  Stress  Intensity  Factor 

As  mentioned  at  the  outset,  the  BJS-based  dynamic  stress  intensity  factor 
Kj  did  not  achieve  a  range  of  values  deemed  critical  for  fracture  until  well 
after  the  experimentally-determined  (Nl)  fracture  initiation  time.  The  ob¬ 
servations  made  above  about  Fig.  8  show,  of  course,  why  this  was  so,  and  we  now 
turn  to  the  FE  results  to  estimate  a  critical  value  of  Kj.  Extrapolation  of 

the  Kj-curve  in  Fig.  7  out  to  the  NL-estimated  fracture  initiation  time 

3/2 

t=70  psec  gives  Kj-144  MN/m  .  This  value  is  higher  than  the  critical  range 
adopted  in  BJS.  However,  that  range  was  based  on  zero- veloci ty  extrapolations 
of  intensity  factor  data  found  for  growing  cracks,  e.g.  (Shockley,  Kalthoff, 
Klemm  and  Winkler,  1983).  Similar  experiments  which  examine  near  fracture 
initiation  have,  however,  been  performed  for  4340  steel  (Zehnder  and  Rosakis, 
1986).  There,  Kj  is  found  to  be  in  the  range  200-250  MN/nf^  at  5  usee  after 
fracture,  but  to  drop  rapidly  afterwards. 

Brief  Discussion 

This  paper  presented  another  step  in  the  full  transient  study  of  a  dynamic 
tear  (DT),  or  three-point  bend,  brittle  fracture  test  modeled  as  a  process  of 
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plane  strain.  As  in  a  previous  seep  (3rock,  Jolles  and  Schroedl .  1985 > .  a 
major  goal  of  chis  paper  was  to  accurately  determine  the  role  of  elastic  wave 
arrivals  at  the  test  specimen  notch  in  generating  the  stress  state  necessary 
for  dynamic  fracture  to  occur . 

An  FE  code,  ABAQUS ,  with  capabilities  for  treating  exact  dynamic  equations 
was  first  compared  with  the  analytical  solution  which  is  exact  for  times  prior 
to  the  arrival  of  wave  reflections  The  FE  results  were  then  used  to  examine  rn 
test  model  behavior  during  the  interval  in  which  the  specimen  notch  receives 
first  reflections  from  the  notched  side  of  the  specimen.  This  behavior 
established  clearly  that  it  is  such  reflections,  and  not  the  initial  impact 
waves,  which  generate  the  stress  state  critical  for  dynamic  fracture. 

The  time  history  of  the  dynamic  stress  intensity  factor  could  also  be 
obtained  graphically  from  the  FE  results,  and  extrapolation  of  this  history 
to  the  experimentally-determined  (Nash  and  _ange ,  1969)  fracture  initiation 
time  yielded  an  estimate  for  the  critical  value.  This  value  was  in  qualitative 
agreement  with  values  obtained  from  ocher,  similar,  tests. 

With  the  importance  of  the  first  reflected  waves  thus  established,  future 
work  will  anlayze  the  test  model  results  near  the  actual  time  of  fracture 
initiation.  More  importantly,  perhaps,  future  work  will  incorporate  the 
possibility  of  material  nonlinearity,  so  that  the  test  model  interpretation 
no  longer  need  be  restricted  to  purely  brittle  fracture. 
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Figure  2  Reflection-free  solution  geometry 
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otch(crack  tip)  at  9.80x10  b  sec 
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FINITE  ELEMENT  RESULTS 
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